Metastable tight knots in a worm-like polymer 
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Based on an estimate of the knot entropy of a worm-like chain we predict that the interplay of 
bending energy and confinement entropy will result in a compact metastable configuration of the 
knot that will diffuse, without spreading, along the contour of the semi-flexible polymer until it 
reaches one of the chain ends. Our estimate of the size of the knot as a function of its topological 
invariant (ideal aspect ratio) agrees with recent experimental results of knotted dsDNA. Further 
experimental tests of our ideas are proposed. 



While everyday experience suggests that knots are very 
common in long linear strings of any kind, and they get 
quite tight whenever a string is not carefully handled, the 
study of knots in polymers concentrated almost exclu- 
sively on closed loops. It is understandable in the sense 
that knots are not mathematically well defined for open 
strings. However, if the string is long enough, while the 
knot occupies a short fragment of it far from the ends, 
then distinguishing between different knots, or between 
knots and no knots, becomes sufficiently unambiguous. 

The recent achievement in the theory of knotted loops 
[1-3] is the idea of knot localization. As simulations show 
[3] , when a polymer loop with a knot is placed in a good 
or 0-solvent, it typically adopts a conformation in which 
most of the polymer forms a long unknotted loop, while 
the knot gets somewhat tightened in a small part of the 
contour. A theoretical explanation of the knot localiza- 
tion phenomenon is given [2] for the so-called flat knots, 
confined in a thin slit between two planes. In this case, 
one can consider a two-dimensional network which cor- 
responds to any configuration of the knot by identifying 
chain crossings as the cross-links. This mapping allows 
one to resort to the sophisticated theory of 2D networks 
with excluded volume [4, 5] and conclude that the most 
entropically favorable network is obtained when one of 
the knot arcs is made long at the expense of all others, 
thus localizing the knot. 

In the present note, we would like to use much sim- 
pler hand- waving arguments to show that, at least for 
a worm-like polymer, there exists a local (metastable) 
minimum of free energy corresponding to a tight state of 
the knot even when the chain itself is open (not a loop). 
That means that if we intentionally tie a sufficiently tight 
knot somewhere on a very long polymer chain, the knot 
will spontaneously shrink or expand to a well-defined size 
and then, on much longer time scales, it will diffuse along 
the polymer (by polymer self-reptating through the knot) 
until, finally, the knot is released through the chain end. 
We expect that this knot will diffuse along the polymer 
as a soliton, in the sense that its size will remain rela- 
tively stable as it diffuses over large distances. The size 
of such a solitary knot depends on the complexity of the 



knot and on the persistence length of the polymer; we 
estimate that the knot will tighten to a size smaller than 
the persistence length of the polymer. 

To put forward our argument, we employ the knot en- 
tropy estimates based on the idea [6], similar to that in 
reptation theory [7] , that non-crossing constraints im- 
posed on the chain in the knot can be self-consistently 
described by confining the chain in an effective tube. 
Physically most natural construction of such a tube cor- 
responds to making maximally inflated or, equivalently, 
shortest length tube consistent with the given knot topol- 
ogy. This approach to estimate knot entropy was in- 
vented by us [6]. It was also used in a number of other 
contexts and widely popularized under the name "ideal 
knots" [8]. The configuration and aspect ratio, p, of the 
"ideal" tube represent topological invariants of the knot. 
This is sketched in the figure 1. 

Within the framework of this approach, we imagine 
that a tight knot, characterized by ideal aspect ratio p, 
has been tied in a very long polymer, and that currently 
the degree of tightening of the knot is such that the di- 
ameter of its self-consistently confining tube is D. The 
length of the tube is pD and the size of the knot in space 
is R, such that pD x D 2 ~ R 3 or 
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Throughout this paper, we drop all numerical coeffi- 
cients, emphasizing only the scaling aspects of our con- 
siderations. We want to estimate the free energy of the 
knot as a function of D or R. 

Although the original estimates of knot entropy were 
designed for a Gaussian chain [6], they can be readily 
adapted to a worm-like polymer. Consider a knot with 
D < £, where £ is the persistence length of the polymer. 
In this case, the polymer is quite tightly confined in the 
tube and can only wiggle a little bit around its centerline. 
Therefore, the length of polymer within the knot is close 
to pD, and its free energy consists of bending energy and 
confinement entropy. We estimate the bending energy as 
TpDij R 2 , assuming temperature T is expressed in en- 
ergy units (ks = 1), and assuming that the radius of 
curvature of the tube is about the knot size R, which is 
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FIG. 1: (Color online) A knot is tied on a polymer chain. 
Knot size in space is R, chain within the knot is confined to 
a self-consistent tube of diameter D. Length of this tube is 
pD, where p is the topological invariant of the knot. 



natural for an "ideal" knot. We further assume that the 
confinement entropy is the same as that for a straight 
tube for which case it was computed by Odijk [9] and 
turns out to be about unity per every so-called deflection 
length A ~ (-D 2 ^) 1 / 3 ; this results in the free energy con- 
tribution about TpD/X ~ Tp(Dfe) 1 ^. Thus, using also 
(1), we obtain the following free energy estimate: 
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where AF is the free energy penalty for forming the knot 
(taking the reference free energy to be that of the unknot- 
ted chain). Obviously, it has a minimum at 
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Notice that the resulting optimal D* meets the condition 
D* < I, so our estimate is self-consistent in this respect. 

Of course, our result applies only as long as D* > d, 
where d is the thickness of the polymer itself. Further- 
more, in the practically important case of a charged poly- 
mer, such as DNA, we should also require that tube diam- 
eter exceeds Debye screening length, D* > r s . In general, 
we can write roughly p < (£/(d + r s )) 2 . More complex 
knots get so tight that their further collapse is stopped 
by either the excluded volume or electrostatic repulsion. 
As an example, for the dsDNA, the ratio £/(d + r s ) w 20 
under physiological conditions [10], while all knots with 
7 or fewer crossings on the projection have p about 30 
or less [11]. Therefore, in practice the condition on knot 
complexity p < (£/(d + r s )) 2 is not very restrictive. 

Let us now try to understand the physical meaning of 
our result (3), because at the first glance it might seem 
counterintuitive. Indeed, the optimal D* results from the 
competition of two factors, each of which, as it seems, 
disfavors tightening! One, chain bending energy, obvi- 
ously favors more loose states of the knot, or increasing 
D. The other one, however, related to the confinement 
entropy, favors tube widening if the chain length in the 
tube is fixed. In our case, when a knot tightens, it re- 
duces the length of the chain, pD, confined in the knot. 
In other words, the part of the chain that remains in the 



tube suffers more when D decreases, but the other part 
of the chain gets completely free of restriction, and that 
factor wins the whole game. That means, what really 
tightens the knot is the entropy gain of the chain tails 
outside of the knot. Alternatively, we can say that the 
knot gets compressed by the pressure of Rouse modes 
(or bending phonons) of the outside chain tails. Similar 
force can be observed during the translocation of a long 
polymer through a narrow channel across a membrane; 
while two long polymer ends are outside the membrane, 
their entropic favorability results in a force stretching 
the polymer portion inside the channel, and, accordingly, 
compressing the membrane. Notice that while in this pla- 
nar geometry this force increases logarithmically with the 
length of the polymer tails, in our case of two long poly- 
mer ends sticking out of a compact knot, the force on the 
knot is independent of the length of these ends. 

How tight should be the knot in the first place in or- 
der for our mechanism to take over and to bring the knot 
to its metastable size R*l In other words, how wide 
is the basin of attraction of our metastable free energy 
minimum? We argue that the knot should be initially 
tightened to the state in which tube diameter is smaller 
or about chain persistence length £. Indeed, when knot 
tube is wider, D > £, the chain inside the tube is roughly 
Gaussian, with blobs of size D. Each blob contains about 
D 2 /£ of polymer contour length, and the entire tube con- 
tains polymer length pD 2 /£. At the same time, confine- 
ment entropy is about unity per blob, which results in 
overall confinement entropy of about p, independent of 
D. Thus, our entropic knot tightening effect does not 
work if the tube is wider than persistence length, and 
it comes into play only when the knot is prepared in a 
compact enough state such that D < £. To achieve this, 
the knot should be initially prepared by the pulling the 
ends with force / > /*, where /* ~ T/l. 

Let us discuss now the dynamics aspect of the situa- 
tion. Imagine once again that a sufficiently tight knot 
was initially tied in the polymer, similar to how it was 
done in the experimental work [12] (see also [13]) with 
DNA. Suppose now that the chain is released and is free 
to move. We predict then that the knot tightens spon- 
taneously within a time which is roughly independent of 
the total chain length L. After that, the knot will diffuse 
along the chain in pretty much the same way as it was 
observed experimentally [12, 13] and numerically [14] for 
the stretched chain (even though in the experiment of ref. 
[12] the chain ends were held at fixed separation through- 
out the diffusion process). As regards the knot diffusion 
coefficient along the chain, it was shown in the work [12] 
that it can be quite accurately expressed in terms of the 
friction coefficient, £, of a polymer with length pD and 
diameter d, moving in the viscous solvent inside a tube 
of diameter D: £ = \ n ^ r j^ pD. The diffusion coefficient 
is then determined by Einstein's relation as T/£, and the 
full relaxation time of knot diffusion to the chain end is 
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FIG. 2: (Color online) Knot size 7?* against the stretching 
force /. Force is indicated in piconewtons and knot size in 
nanometers, both under the assumptions that I = 50 nm, and 
that there are no any numerical coefficients to worry about. 



about L 2 (/T. In these formulae, r\ is assumed to be the 
bulk viscosity of the solvent. 

The above description of a soliton-like knot diffusing 
as a whole along the chain holds as long as diffusion time 
is shorter than the time of thermally-activated loosening 
of the knot. Indeed, we argued that the knot gets to its 
metastable size R* only when it is compact enough to 
begin with, such that D < £ and, therefore, a free energy 
barrier exists at D w £ (knots with D > £ will loosen up 
spontaneously). Using formula (2), we can estimate the 
barrier height as AF\ D/e=1 - AF\ D=D , ~ T (p-p 5 / 6 ). 
As one could have expected, this barrier is very high for 
complex knots, but even for p « 30, one gets a bar- 
rier of about 10T which should be sufficient to keep a 
knot locked while diffusing over a distance of few hun- 
dred persistence lengths, typical of DNA manipulation 
experiments. For smaller values of p, the barrier may 
be too small to stabilize the tight knot and without the 
applied stretching it will spread due to thermal fluctua- 
tions. 

As a corollary to our result, let us mention the follow- 
ing rather unexpected prediction. Let us take the poly- 
mer with the knot tightened according to our mechanism, 
and now let us gently pull the chain ends by a weak force. 
The applied stretching will suppress transverse fluctua- 
tions of the open chain ends. But these fluctuations are 
exactly the reason why the knot was tightened in the first 
place. Therefore, when these fluctuations are suppressed 
by the applied force, the knot swells - instead of fur- 
ther tightening which one could have naively expected. 
Of course, at the larger forces the loosening of the knot 
stops and normal tightening takes over. 

Let us support this physical argument by a little calcu- 
lation. If the chain of contour length L is stretched by a 
weak force /, it represents the succession of Pincus blobs, 
each involving the number of persistence lengths g such 
that fig 1 / 2 ~ T; the stretching free energy is about T 
per blob, or TL/gi ~ L£f 2 /T (see, e.g., [15] for further 
details). In our case, only the chain outside of the knot is 



subject to this stretching effect, which yields the overall 
free energy 

subject to optimization with respect to D (see Eq. (6) 
below); here, the term oc L is large, but can be dropped 
as independent of D. Formula (4) remains valid as long 
as Pincus blob is larger than persistence length ft/T < 1. 

If we stretch the outside chain even further, beyond the 
Gaussian regime, it crosses-over to the so-called Marco- 
Siggia or Odijk regime in which stretching free energy is 
about T per Odijk deflection length. In this regime, 
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once again, this has to be optimized with respect to D. 

Notice that in both equations (4) and (5) we have ne- 
glected the effect of force on the chain part inside of the 
knot. This is justified as long as the amount of trans- 
verse fluctuations of the outside chain, D out (f) 1 which is 
either the blob size in Pincus regime or the tube diame- 
ter in Odijk regime, remains smaller than the optimized 
tube diameter inside the knot, D* . At still larger forces, 
the chain inside the knot is just as stretched as outside. 
Therefore, we can summarize all results of optimization 
in the following way: 
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These results are sketched in figure 2. In accordance with 
our qualitative argument, the application of weak force 
loosens the tight knot instead of tightening it further! 

In practice, for the case of dsDNA, since i » 50 nm 
and T ~ 4 pN x rtm, and taking p w 30, we get that 
both cross-overs fl/T ~ 1 and f£/T ~ p 1 / 3 are within 
an experimentally feasible force range (about a few tenth 
and about a few piconewtons, respectively). 

Speaking about experimental tests of our theory, we 
should of course consider first the experiment by Quake 
et al [12, 13] in which the knots were tied on DNA by the 
use of optical tweezers and then the knot motion along 
DNA was observed. Our theory indicates that the op- 
timal tube diameter in the knot, D*, does depend in a 
certain way on the knot complexity, on p. In the exper- 
iment, authors were unable to observe directly the knot 
and measure its size R, but they were able to detect the 
amount of fluorescence coming from the knot region in ex- 
cess of the fluorescence from linear DNA. That means, in 
our notations, that authors measured quantity pD* — R* , 
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FIG. 3: (Color online) The excess DNA length in the knot, 
pD* — R* against p — p 1 ^ 3 . The meaning of these variables is 
explained in the text. Data points are taken from the work 
[12, 13]. Linear fit corresponds to the assumption that D* 
does not depend on the knot type, on p. Curved line corre- 
sponds to our theory. 

because pD* is the length of DNA within the knot (pro- 
portional to the amount of fluorescence), while R is the 
length of DNA that would have been in that place if 
there were no knot. Given the normalization condition 
(1), one could write pD* - R* = (p - p 1 / 3 ) D* . Authors 
of the work [12], assuming a'priori that D is indepen- 
dent of p, plotted pD* — R* against p — p 1 ^ 3 , fitted the 
data to the linear function, and interpreted its slope as 
the tube diameter D* . By contrast, our theory predicts 
that D* does depend on knot complexity p, such that 
pD* -R* = {p 1 ' 2 - p- 1 ^) I. This is illustrated in the 
figure 3, where both the linear fit and our theory results 
are plotted along with the data points reproduced from 
the work [12]. At present, it is impossible to decide which 
theory is a better fit to the data. 

To further test our predictions, one should prepare the 
knot and then consider what happens to it if the force is 
subsequently switched off or significantly reduced. An- 
other test of our theory would be to see how the knot 
diffusion time, determined by the friction coefficient £, 
depends on the knot type through D. Yet another possi- 
bility is to use granular macroscopic chain experiments, 
as described in [16]. In the movie available on the site 

http : //cnls . lanl . gov/~ebn/research/pics/chain .mov 

the knot does move along the chain like a soliton, it does 
not increase in size while moving, in accord with our the- 
ory. It would be interesting to make a detailed statistics 
on this subject, although we should also put in question 
the applicability of our arguments based on worm-like 
chain model to the granular chain examined in [16]. 

To conclude, we have shown that a sufficiently complex 
knot in a worm-like polymer will attain a well-defined 
compact size which is smaller that the persistence length. 
This configuration is metastable and the knot will diffuse 
along the chain as a soliton, whether the ends of the chain 
are subjected to external force or are free to move. No- 



tice that since the results depend crucially on the elastic 
properties of the polymer on length scales comparable to 
and smaller than the persistence length, the effects de- 
scribed in this work are quite sensitive to the choice of 
the polymer model and can not be captured by locally 
stiff (e.g., freely jointed or freely rotating) chain models. 
Likewise, solitary knots will not form in Gaussian chains 
for which the persistence length £ and the thickness of 
the chain d coincide. Such knots are expected to form in 
semi-flexible polymers such as double stranded DNA for 
which (a) separation of length scales exists, L ^> £ 3> d 
and (b) elastic behavior on length scales smaller than £ 
was demonstrated in single molecule experiments [17, 18]. 
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